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^ I Abstract. 

^ ■ The power spectrum of transmitted flux in the Lya forest in spectra of distant 

I quasars gives us an information about matter power spectrum on the smallest 

■ spatial scales which is very important for testing of different models of dark matter. 
I Only a few independent samples of spectra obtained on diferent instruments have 

^ been used for this purpose, thus using of additional independent sample is useful 

I for verification an existing results. We used the data from the 2dF (Two-degree 

r \ • Field) QSO Redshift Survey (2QZ) to obtain the mean transmitted flux for three 

I different redshift bins. After visual inspection and rejection of spectra with broad 

■ absorption lines, damped Lya systems and low signal-to-noise ratio spectra our 
'p^l final sample contains 655 quasars with redshift range 2.3 < z < 2.9. We present 
O ' the results on composite spectra, determination of continuum level for different 
^ I redshift ranges and the mean transmission in the Lya forest as a function of redshift 
^ ■ within the range 2.1 < Zfo^est < 2.5. The obtained redshift-dependence of the mean 

transmission agrees well with the results of other authors. These results can be 
used for calculation of the flux power spectrum of Lya forest. 

> ■ 

o : 

\Q , Introduction 

m ■ 

^Sj ■ The prominent absorption features of the spectra shortward of the Lya emission line are 

\^ . called Lya forest and were first discovered by Gunn and Peterson in 1965 like a small decrement 

I in the spectra [Gunn et al, 1965]. The high resolution measurements have been largely combined 

■ with hydrodynamical simulations and semi-analytical models. The results of these simulations 
• • . have shown that the neutral gas responsible for the absorption is in a relatively low density, 

I smooth environment, which implies a simple connection between the gas and the underlying 

^ ' dark matter. Thereby the fluctuations in the Lya optical depth should reflect the density of 

■ the dark matter distribution on scales larger than a "filtering" scale related to the Jeans length. 
This paradigm for the origin of the Lya forest has led to considerable interest in using quasar 
absorption spectra to study the dark matter distribution. Of particular interest is thereby the 
possibility to probe the density fluctuations of the matter with the flux power spectrum (FPS) 
of quasar absorption lines. 

Only a few independent samples of spectra obtained on diferent instruments have been 
used for this purpose. In [Kim et al, 2004] LUQAS sample consisting of 27 spectra obtained 
with the ESO Ultra- Violet Echelle Spectrograph (UVES) on VLT, Paranal, Chile, was used 
for calculations of the Lya forest FPS for three redshift bins with (z) = 1.87, {z) = 2.18 
and (z) = 2.58. This sample has spectral resolution ~ 45 000, signal-to-noise ratio larger 
than 25 and wavelength coverage 3000 - 11000 A. In [Viel et al, 2004] the same sample and 
additional data from [Croft et al, 2002] together with a large suite of high-resolution large 
box-size hydrodynamical simulations were used to estimate the mean transmitted flux and 
FPS for three redshift bins with {z) = 2.125, {z) = 2.44 and {z) = 2.72. In [Cappetta et al, 
2010] the sample of 15 quasar spectra, observed with UVES was used to obtain the auto- and 
cross-correlation functions of the Lya forest with median redshift {z) ~ 1.8. This sample has 
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signal-to- noise ratio ~ 3 - 12. On the other hand, McDonald et al, [2006] used 3035 quasar 
spectra from the Sloan Digital Sky Survey (SDSS), obtained with multi-fibre spectrograph that 
has resolution ^ 2000, wavelength coverage 3800 - 9200 A and is placed at Apache Point 
Observatory, New Mexico, USA. The authors presented FPS of the Lya forest for 9 redshift 
bins within the range of 2.2 < z < 3.8. SDSS spectra have much lower resolution and signal-to- 
noise ratio than UVES ones, but they cover much larger redshift range and count much more 
spectra that plays a significant role in the statistical errors reduction. 

Thereby, it is useful to involve the independent data samples obtained on different equip- 
ment in the large-scale structure investigations with Lya forest. 

One of the main steps in procedure of Lya forest FPS and correlation function calculations 
is estimation of the mean transmitted flux. In the present work we use the independent data 
from the 2dF (Two-degree Field) QSO Redshift Survey (2QZ) to obtain the mean transmitted 
flux for three different redshift ranges. 

The sample 

The 2dF facility used multi-fibre spectrograph on the 3.9-m Anglo-Australian Telescope, 
which is capable of observing 400 objects simultaneously over a 2° diameter field of view. 
Two separate but linked large redshift surveys were the initial main projects with 2dF: one 
for 250,000 galaxies and one for 30,000 color-selected quasars covering redshifts up to 4. Two 
spectographs provide spectra with resolutions of between 500 and 2000, over wavelength range 
of 4400 - 11000 A. The main survey regions are two declination strips, one in the southern 
Galactic hemisphere spanning 80°xl5° around the South Galactic Pole, and the other in the 
northern Galactic hemisphere spanning 75° x 10° along the celestial equator; in addition, there 
are 99 fields spread over the southern Galactic cap. The survey covers 2000 deg^ [Lewis et al., 
2002]. 

Our initial sample consists of 1,961 non fiux-calibrated spectra and covers the redshift range 
of 2.3 < 2; < 3.3. For calibration of the spectra we used the 2dF relative efficiency curve from 
[Lewis et al, 2002]. Calibration included the average telluric lines and fibre absorption. Spectra 
with low signal-to-noise ratio, non-quasar spectra, spectra with broad absorption line features 
(BAL quasars) and the damped Lya systems (DLAs) were removed from the initial sample. The 
DLAs are the systems with a number density of >2x 10^° cm~^ that is three orders of magnitude 
higher than for the Lya forest systems [Meiksin et al, 2009] and appears in the spectrum as 
a wide and deep absorption feature in the Lya forest region. Thus such spectra with DLAs 
can affect the value of the mean transmission. BALs from material along the line of sight to 
the quasar are generated by high-velocity, structured outflows in the central regions of quasars. 
Absorption signatures of these outflows are, by definition, at least 2000 km s~^ wide and are 
blueshifted [Hall et al, 2002], [Gibson et al, 2010]. The spectra containing BALs can heavily 
obscure the continuum in the resulting composite spectum. After calibration and cleaning the 
resulting sample contains 655 spectra with redshifts within the range of2.3<z<2.9. 

The Method 

Generating the composite spectra 

Subtle global spectral properties and features can be studied by combining large numbers 
of quasar spectra into composites. Such composite spectra have a high signal-to-noise ratio and 
reveal weak emission features that are rarely detectable in individual spectra. This method is 
used by many authors for both high and low signal-to- noice data sets (see [Telfer et al, 1997], 
[Brotherton et al., 2001], [Scott et al, 2004], [Francis et al, 1991] and others). We have made 
composite spectra to obtain the value of the mean transmitted flux for each of them. 

Before making composites all spectra were smoothed with 3-point unweighted sliding- 
average smooth, rebinned to the rest frame and flux normalized. Normalization procedure 
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involved the next steps: estimation of the mean flux within the range of 1450 A < Xrest < 1470 A 
for each spectrum and division each spectrum by this value. This spectral range was chosen as 
the most free one from emission lines. 

We have compiled three composite spectra as arithmetic means for redshift ranges of 2.3 < 
z < 2.5, 2.5 < 2; < 2.7 and 2.7 < 2; < 2.9. The resulting composite spectra cover the rest 
frame wavelength range of 1056 - 1855 A. It is worth to note that the arithmetic mean of the 
power-law spectra will not, in general, result in a power law with a mean exponent but Vanden 
Berk et al, [2001] has shown that the difference in the spectral exponent for mean arithmetic 
and mean geometric composite spectra is insignificant. Such a difference is not important for 
our purpose. 

Continuum fitting 

The estimation of the continuum level is a very delicate step in the process of spectra 
reduction. There is no specific and adopted method for continuum fitting procedure in the 
literature, although different methods had been used (see [Press et al, 1993], [Cappetta et al., 
1993], [Vanden Berk et al, 2001] and others). The problem of fitting the quasar continuum is 
complicated by the fact that there are essentially no absorption-line-free regions in the spectrum. 
We have determined the continuum level for whole composite spectrum using a set of regions 
which are free from emission features (judged by eye). The regions that satisfy this are 1281 - 
1291 A, 1321 - 1327 A, 1348 - 1361 A and 1438 - 1468 A. Continuum level was fitted with a 
power-law / = ^A". The composite spectrum for 2.7 < z < 2.9 with determined continuum is 
shown in Figure 1. The uncertainties for the normalizing coefficient A and for the slope a can 
affect the results reported in this paper on the mean transmitted flux function and should be 
considered as a source of possible systematic errors in addition to quoted statistical errors of 
the results. 
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Figure 1. The composite spectrum for the 2.7 < 2; < 2.9 range (solid line with data points) 
and estimated continuum (dashed line). 

The mean transmitted flux 

The matter density fluctuations can be calculated from the FPS of Lya forest which in 

turn is obtained using the mean transmitted flux in this region of spectrum. The flux power 
spectum is the measure of the variance in the amplitude of the Fourier transform coefficients of 
the transmitted flux. 

The transmitted flux of a quasar spectrum is given by F'{X) = Fccmti^)e~'^^^\ where 
-Pcont(A) is the continuum, e"'^ is the transmission and r(A) is the opacity or optical depth. For 
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the spectrum reduced by the continuum fitting: 

F(A) = e-"W+n(A), (1) 

where n(A) describes the stochastic noise. Averaging ([T]) over the density fluctuations (over 
all pixels in a redshift bin) one can obtain ^(A) = (e"'^^^^). Generally speaking, this value is 
position dependent, as it depends on the photoionization rate of HI and the temperature of 
the intergalactic medium. If one assumes that -F(A) is a constant, this value then is the mean 
transmitted flux F{z) which is a function of the redshift. 

The mean transmitted flux can be determined through the folowing considerations. The 
flux density in the i-th pixel in a single spectrum within the forest range is 

//(A.e.t, z) = A^C{Xrest)[F{z) + SF^] + n|, (2) 

where j stands for spectrum, is the normalizing coefficient, C{Xrest) is the continuum, F(z) 
is the mean transmitted flux, SF- is the fluctuations of the mean transmission and is the 
noise. Averaging over a small redshift bin (a set of spectra) for the composite spectra we obtain 

fi = C{\rest)F{z), (3) 

as {6Fi) = over all lines of sight and (n-j) = over all spectra. We determine the continuum 
level inside the forest region (1056 - 1202 A) as 



C = A {XrestT + ci exp 



(Xrest - 1066.66) 



21 



2cl 



+ C3 exp 



{Kest - 1123.00) 



21 



2cl 



+ 



+ C5 exp 



{Xrest - 1175.70)^ 



2ci 



+ C7 exp 



{Xrest - 1215.67)^ 

2cl 



, (4) 



where A and a are the mean continuum coefficients, the wavelenghts of 1066.66 A, 1123.00 A, 
1175.70 A and 1215.67 A are the wavelengths of Arl, FellL, CIII* and Lya spectral lines 
respectively. 



Results 

For three redshift bins 2.3 < z < 2.5, 2.5 < z < 2.7 and 2.7 < z < 2.9 we derive continuum 
slopes of a ~-0.6, a ~~0.5 and a ~-l.l respectively. The resulting mean transmitted flux and 
parameters of the lines are presented in Table 1. Estimated la error for F includes statistical 
errors of the fitted data. For calculation of the mean transmitted flux mean continuum param- 
eters were adopted. For illustration in Figure 2 the Lya forest region with the obtained fitting 
curve are shown for 2.5 - 2.7 redshift bin. 

Table 1. Line parameters used for continuum fitting and the mean transmitted 
flux for three redshift bins. 



Quasars redshift 
bin 


^ falsest 


Arl 


Felll: 


CIII* 


Lya 


F 


cl 


c2 


c3 


c4 


c5 


c6 


c7 


c8 


2.3-2.5 


2.13 ± 0.29 






0.2 


30.5 


0.1 


11.3 


1.5 


15.2 


0.83 ± 0.04 


2.5-2.7 


2.29 ± 0.30 


0.3 


22.6 


0.04 


6.8 


0.2 


9.2 


1.7 


16.2 


0.78 ± 0.03 


2.7-2.9 


2.47 ± 0.31 


0.2 


10.5 


0.1 


6.4 


0.2 


7.0 


1.5 


16.2 


0.69 ± 0.03 



In Figure 3 we present the comparison of the obtained mean transmitted flux with the 
results of other authors. Results of McDonald et al, [2000], where eight high-resolution and 
high signal-to-noise ratio spectra obtained with HIRES instrument were used, are shown with 
dots. Results of Viel et al, [2004] where the mean transmitted flux was obtained for three 
redshift bins using the spectra from the LUQAS sample, are shown with triangles. Two dashed 
lines show the upper and lower limits for 7 in the power-law describing the mean transmitted flux 
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Figure 2. Lya forest region (data points) with spectral lines (pointed with arrows) and obtained 
fitting curve (solid line) for composite spectrum of the 2.5 < z < 2.7 redshift bin. 
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Figure 3. Obtained mean transmitted flux compared with other authors results (for the 
explanation see the text). 

evolution described by Press et al., [1993] as function f (z) = A[l + z)^~^^ , where the coefficients 
7 = 2.46 ± 0.37 and A = 0.0175 - O.OO567 ± 0.0002 were obtained using this power-law to 
fit the data. Authors used 29 spectra from the Schneider-Schmidt-Gunn (SSG) sample. Our 
results are shown with squares. The error bars in z direction show the Lya forest range in each 
bin. 

Conclusions 

Independent result of calculation of the mean transmitted flux was obtained for three 
redshift bins using the low S/N and low resolution data from the 2dF QSO survey. Mean trans- 
mitted flux is F{z) = 0.83 ± 0.04, F{z) = 0.78 ± 0.03, F{z) = 0.69 ± 0.03 for the redshift bins 
of the Lya forest Zforest = 2.13, 2.29 and 2.47 respectively. The obtained redshift-dependence 
of the mean transmission agrees well with the results of other authors (e.g. McDonald et al, 
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[2000], Viel et al, [2004], Press et ai, [1993]). It is notable that the derived results on the 
mean transmission have lower values in comparison with other authors. This can be explained 
by the uncertainties of the continuum level estimation and usage of the simple power-law for 
this purpose. Besides, we have shown that the data from 2QZ is reasonable for the Lya forest 
investigations. These results can be used for calculation of the flux power spectrum of Lya 
forest. 
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